Domain 3a mutation of VPS33A suppresses larval arrest phenotype in the loss of VPS45 in Caenorhabditis elegans.

The Sec1/Munc18 (SM) protein VPS45 is a key regulator of SNARE-mediated membrane fusion in endosomal trafficking, but its precise role remains unknown. To understand the function of VPS45 in vivo , we performed a genetic suppressor screen in Caenorhabditis elegans . We found that the temperature-sensitive lethality caused by the loss of VPS-45 can be suppressed by a mutation in another SM protein, VPS33A. The VPS33A M376I mutation is located in domain 3a, which is predicted to be essential for SNARE complex assembly. These results highlight the functional importance of domain 3a in endosomal SM proteins and its role in specific membrane fusion.

Metagenomic analyses of 7000 to 5500 years old coprolites excavated from the Torihama shell-mound site in the Japanese archipelago.

Coprolites contain various kinds of ancient DNAs derived from gut micro-organisms, viruses, and foods, which can help to determine the gut environment of ancient peoples. Their genomic information should be helpful in elucidating the interaction between hosts and microbes for thousands of years, as well as characterizing the dietary behaviors of ancient people. We performed shotgun metagenomic sequencing on four coprolites excavated from the Torihama shell-mound site in the Japanese archipelago. The coprolites were found in the layers of the Early Jomon period, corresponding stratigraphically to 7000 to 5500 years ago. After shotgun sequencing, we found that a significant number of reads showed homology with known gut microbe, viruses, and food genomes typically found in the feces of modern humans. We detected reads derived from several types of phages and their host bacteria simultaneously, suggesting the coexistence of viruses and their hosts. The food genomes provide biological evidence for the dietary behavior of the Jomon people, consistent with previous archaeological findings. These results indicate that ancient genomic analysis of coprolites is useful for understanding the gut environment and lifestyle of ancient peoples.

Bioarchaeological study of ancient Teotihuacans based on complete mitochondrial genome sequences and diet isotopes.

BACKGROUND: The Teotihuacan civilisation was the largest one in ancient Mesoamerica. The Teotihuacan city was born in the north-eastern Basin of Mexico around the second century BC, reached its peak in the fourth century AD, and had cultural influence throughout Mesoamerica. At its peak, the size of the city reached more than 20 km2, and the total population is estimated to have increased from 100,000 to 200,000. However, knowledge of the genetic background of the Teotihuacan people is still limited. AIM: We aimed to determine the mitogenome sequences of the Teotihuacan human remains and compare the ancient and present Mesoamericans. In addition, we aimed to identify the food habits of ancient Teotihuacans. SUBJECTS AND METHODS: We determined the mitogenome sequences of human remains dated to 250-636 cal AD using target enrichment-coupled next generation sequencing. We also performed stable isotope analysis. RESULTS: We successfully obtained nearly full-length sequences newly unearthed from a civilian dwelling in the Teotihuacan site. Teotihuacan mitochondrial DNA was classified into the haplogroups in present and ancient Mesoamericans. In addition, Teotihuacan individuals had a diet dependent on C4 plants such as maize. CONCLUSION: Genetic diversity varied among the Teotihuacans.

A calf with hind limb paralysis and dysstasia and a genome sequence analysis of an isolated Clostridium perfringens toxinotype E strain.

Clostridium perfringens toxinotype E infections are rare in calves, and the development of intestinal lesions were commonly observed. In 2012, a 6-day-old calf in Japan exhibited swelling with emphysema on the right gluteal region, sudden paralysis of the hind limb and dysstasia. A pathological examination revealed myositis of the gluteal muscle and neuritis of the ischiatic nerve. C. perfringens type E strain CP118 was isolated from the affected muscle. However, the intestinal symptoms and lesions that commonly develop in type E infections in calves were not detected in the present case. Genome analyses revealed that CP118 possessed 16 virulence-related genes, including enterotoxin, and was closely related to other type E and F strains. Particularly, CP118 was more closely related to type E strains from humans, including a food poisoning case, than calf isolates, suggesting its potential to cause food poisoning in humans and, thus, its importance as a potential risk to public health. Since CP118 did not possess the reported toxin genes associated with neuropathy, pyogenic inflammation caused by CP118 and/or other bacteria may have damaged the ischiatic nerve, resulting in neuropathy. Alternatively, unidentified CP118 toxins may have caused the neuropathy. This is the first study to report C. perfringens type E infection with peripheral neuropathy. The distribution of all the reported virulence-related genes in the C. perfringens population as well as the details of this rare case will provide further insights into C. perfringens type E infections.

Detection of non-reference porcine endogenous retrovirus loci in the Vietnamese native pig genome.

The Vietnamese native pig (VnP)-a porcine breed with a small body-has proven suitable as a biomedical animal model. Here, we demonstrate that, compared to other breeds, VnPs have fewer copies of porcine endogenous retroviruses (PERVs), which pose a risk for xenotransplantation of pig organs to humans. More specifically, we sought to characterize non-reference PERVs (nrPERVs) that were previously unidentified in the reference genome. To this end, we used whole-genome sequencing data to identify nrPERV loci with long terminal repeat (LTR) sequences in VnPs. RetroSeq was used to estimate nrPERV loci based on the most current porcine reference genome (Sscrofa11.1). LTRs were detected using de novo sequencing read assembly near the loci containing the target site duplication sequences in the inferred regions. A total of 21 non-reference LTR loci were identified and separated into two subtypes based on phylogenetic analysis. Moreover, PERVs within the detected LTR loci were identified, the presence of which was confirmed using conventional PCR and Sanger sequencing. These novel loci represent previously unknown PERVs as they have not been identified in the porcine reference genome. Thus, our RetroSeq method accurately detects novel PERV loci, and can be applied for development of a useful biomedical model.

Prevotella lacticifex sp. nov., isolated from the rumen of cows.

The genus Prevotella plays an important role in polysaccharide degradation and fermentation in the rumen. To further understand the function of the phylogenetically diverse genus Prevotella, it is necessary to explore the individual characteristics at the species level. In this study, Gram-negative anaerobic bacterial strains isolated from the rumen of Holstein cows were identified. Strain R5019T was classified within the genus Prevotella based on 16S rRNA gene sequence-based phylogenetic analysis. The values of 16S rRNA gene sequence similarity, average nucleotide identity and digital DNA-DNA hybridization between strain R5019T and its phylogenetically nearest species Prevotella multisaccharivorax PPPA20T were 89.8, 82.6, and 29.3 %, respectively. The genome size of R5019T was estimated to be ca. 4.19 Mb with a genomic G+C content of 49.5 mol%. The major cellular fatty acids and menaquinones were C15 : 0 anteiso and C17 : 0 anteiso and MK-11 and MK-12, respectively. Succinate, lactate, malate, acetate and formate were produced as the fermentation end products using glucose. Based on phylogenetic, physiological, biochemical and genomic differences between 11 strains and other phylogenetically related Prevotella species, a novel species, Prevotella lacticifex sp. nov., is proposed within the genus Prevotella. The type strain is R5019T (=JCM 34664T=DSM 112675T).

Development of an SSR marker-based genetic linkage map and identification of a QTL associated with flowering time in Eustoma.

Lisianthus (Eustoma grandiflorum) is an important floricultural crop cultivated worldwide. Despite its commercial importance, few DNA markers are available for molecular genetic research. In this study, we constructed a genetic linkage map and to detect quantitative trait loci (QTLs) for important agronomic traits of lisianthus. To develop simple sequence repeat (SSR) markers, we used 454-pyrosequencing technology to obtain genomic shotgun sequences and subsequently identified 8263 putative SSRs. A total of 3990 primer pairs were designed in silico and 1189 unique primer pairs were extracted through a BLAST search. Amplification was successful for more than 1000 primer pairs, and ultimately 278 SSR markers exhibited polymorphism between the two lisianthus accessions evaluated. Based on these markers, a genetic linkage map was constructed using a breeding population derived from crosses between the two accessions, for which flowering time differed (>140 days when grown under 20°C). We detected one QTL associated with flowering time (phenotypic variance, 27%; LOD value, 3.7). The SSR marker located at this QTL may account for variation in flowering time among accessions (i.e., three accessions whose nodes of the first flower were over 30 had late-flowering alleles of this QTL).

Population dynamics in the Japanese Archipelago since the Pleistocene revealed by the complete mitochondrial genome sequences.

The Japanese Archipelago is widely covered with acidic soil made of volcanic ash, an environment which is detrimental to the preservation of ancient biomolecules. More than 10,000 Palaeolithic and Neolithic sites have been discovered nationwide, but few skeletal remains exist and preservation of DNA is poor. Despite these challenging circumstances, we succeeded in obtaining a complete mitogenome (mitochondrial genome) sequence from Palaeolithic human remains. We also obtained those of Neolithic (the hunting-gathering Jomon and the farming Yayoi cultures) remains, and over 2,000 present-day Japanese. The Palaeolithic mitogenome sequence was not found to be a direct ancestor of any of Jomon, Yayoi, and present-day Japanese people. However, it was an ancestral type of haplogroup M, a basal group of the haplogroup M. Therefore, our results indicate continuity in the maternal gene pool from the Palaeolithic to present-day Japanese. We also found that a vast increase of population size happened and has continued since the Yayoi period, characterized with paddy rice farming. It means that the cultural transition, i.e. rice agriculture, had significant impact on the demographic history of Japanese population.

Antimicrobial Resistance Genes in Bacteria Isolated From Japanese Honey, and Their Potential for Conferring Macrolide and Lincosamide Resistance in the American Foulbrood Pathogen Paenibacillus larvae.

American foulbrood (AFB) is the most serious bacterial disease of honey bee brood. Spores of the causative agent Paenibacillus larvae are ingested by bee larvae via brood foods and germinated cells proliferate in the larval midgut. In Japan, a macrolide antibiotic, tylosin, is used as the approved prophylactic for AFB. Although tylosin-resistant P. larvae has yet to be found in Japan, it may emerge in the future through the acquisition of macrolide resistance genes from other bacteria, and bacteria latent in brood foods, such as honey, may serve as a source of resistance genes. In this study, to investigate macrolide resistance genes in honey, we attempted to isolate tylosin-resistant bacteria from 53 Japanese honey samples and obtained 209 isolates from 48 samples in the presence of 1 µg/ml of tylosin. All isolates were Gram-positive spore-forming bacteria mainly belonging to genera Bacillus and Paenibacillus, and 94.3% exhibited lower susceptibility to tylosin than Japanese P. larvae isolates. Genome analysis of 50 representative isolates revealed the presence of putative macrolide resistance genes in the isolates, and some of them were located on mobile genetic elements (MGEs). Among the genes on MGEs, ermC on the putative mobilizable plasmid pJ18TS1mac of Oceanobacillus strain J18TS1 conferred tylosin and lincomycin resistance to P. larvae after introducing the cloned gene using the expression vector. Moreover, pJ18TS1mac was retained in the P. larvae population for a long period even under non-selective conditions. This suggests that bacteria in honey is a source of genes for conferring tylosin resistance to P. larvae; therefore, monitoring of bacteria in honey may be helpful to predict the emergence of tylosin-resistant P. larvae and prevent the selection of resistant strains.

CTCF looping is established during gastrulation in medaka embryos.

Chromatin looping plays an important role in genome regulation. However, because ChIP-seq and loop-resolution Hi-C (DNA-DNA proximity ligation) are extremely challenging in mammalian early embryos, the developmental stage at which cohesin-mediated loops form remains unknown. Here, we study early development in medaka (the Japanese killifish, Oryzias latipes) at 12 time points before, during, and after gastrulation (the onset of cell differentiation) and characterize transcription, protein binding, and genome architecture. We find that gastrulation is associated with drastic changes in genome architecture, including the formation of the first loops between sites bound by the insulator protein CTCF and a large increase in the size of contact domains. In contrast, the binding of the CTCF is fixed throughout embryogenesis. Loops form long after genome-wide transcriptional activation, and long after domain formation seen in mouse embryos. These results suggest that, although loops may play a role in differentiation, they are not required for zygotic transcription. When we repeated our experiments in zebrafish, loops did not emerge until gastrulation, that is, well after zygotic genome activation. We observe that loop positions are highly conserved in synteny blocks of medaka and zebrafish, indicating that the 3D genome architecture has been maintained for >110-200 million years of evolution.

Investigation of the Genetic Diversity of a Rice Core Collection of Japanese Landraces using Whole-Genome Sequencing.

The Rice Core Collection of Japanese Landraces (JRC) consisting of 50 accessions was developed by the genebank at the National Agriculture and Food Research Organization (NARO) in 2008. As a Japanese landrace core collection, the JRC has been used for many research projects, including screening for different phenotypes and allele mining for target genes. To understand the genetic diversity of Japanese Landraces, we performed whole-genome resequencing of these 50 accessions and obtained a total of 2,145,095 single nucleotide polymorphism (SNPs) and 317,832 insertion-deletions (indels) by mapping against the Oryza sativa ssp. japonica Nipponbare genome. A JRC phylogenetic tree based on 1,394 representative SNPs showed that JRC accessions were divided into two major groups and one small group. We used the multiple genome browser, TASUKE+, to examine the haplotypes of flowering genes and detected new mutations in these genes. Finally, we performed genome-wide association studies (GWAS) for agronomical traits using the JRC and another core collection, the World Rice Core Collection (WRC), comprising 69 accessions also provided by the NARO genebank. In leaf blade width, a strong peak close to NAL1, a key gene for the regulation of leaf width, and, in heading date, a peak near HESO1 involved in flowering regulation were observed in GWAS using the JRC. They were also detected in GWAS using the combined JRC + WRC. Thus, JRC and JRC + WRC are suitable populations for GWAS of particular traits.

Draft Genome Sequence of Brucella ceti Isolated in the Western Pacific Ocean.

In 2018, Brucella ceti was isolated from a bottlenose dolphin from the western Pacific Ocean. Here, we report a draft genome sequence of the isolate BD1442 of sequence type 27, which is the only sequence type known to have been isolated from human clinical cases.

Same Locus for Non-shattering Seed Pod in Two Independently Domesticated Legumes, Vigna angularis and Vigna unguiculata.

Loss of pod shattering is one of the most important domestication-related traits in legume crops. The non-shattering phenotypes have been achieved either by disturbed formation of abscission layer between the valves, or by loss of helical tension in sclerenchyma of endocarp, that split open the pods to disperse the seeds. During domestication, azuki bean (Vigna angularis) and yard-long bean (Vigna unguiculata cv-gr. Sesquipedalis) have reduced or lost the sclerenchyma and thus the shattering behavior of seed pods. Here we performed fine-mapping with backcrossed populations and narrowed the candidate genomic region down to 4 kbp in azuki bean and 13 kbp in yard-long bean. Among the genes located in these regions, we found MYB26 genes encoded truncated proteins in azuki bean, yard-long bean, and even cowpea. As such, our findings indicate that independent domestication on the two legumes has selected the same locus for the same traits. We also argue that MYB26 could be a target gene for improving shattering phenotype in other legumes, such as soybean.

TASUKE+: a web-based platform for exploring genome-wide association studies results and large-scale resequencing data.

Recent revolutionary advancements in sequencing technologies have made it possible to obtain mass quantities of genome-scale sequence data in a cost-effective manner and have drastically altered molecular biological studies. To utilize these sequence data, genome-wide association studies (GWASs) have become increasingly important. Hence, there is an urgent need to develop a visualization tool that enables efficient data retrieval, integration of GWAS results with diverse information and rapid public release of such large-scale genotypic and phenotypic data. We developed a web-based genome browser TASUKE+ (https://tasuke.dna.affrc.go.jp/), which is equipped with the following functions: (i) interactive GWAS results visualization with genome resequencing data and annotation information, (ii) PCR primer design, (iii) phylogenetic tree reconstruction and (iv) data sharing via the web. GWAS results can be displayed in parallel with polymorphism data, read depths and annotation information in an interactive and scalable manner. Users can design PCR primers for polymorphic sites of interest. In addition, a molecular phylogenetic tree of any region can be reconstructed so that the overall relationship among the examined genomes can be understood intuitively at a glance. All functions are implemented through user-friendly web-based interfaces so that researchers can easily share data with collaborators in remote places without extensive bioinformatics knowledge.

Fusion of piggyBac-like transposons and herpesviruses occurs frequently in teleosts.

BACKGROUND: Endogenous viral elements play important roles in eukaryotic evolution by giving rise to genetic novelties. Herpesviruses are a large family of DNA viruses, most of which do not have the ability to endogenize into host genomes. Recently, we identified a novel type of endogenous herpesvirus, which we named "Teratorn", from the medaka (Oryzias latipes) genome, in which the herpesvirus is fused with a piggyBac-like DNA transposon, forming a novel mobile element. Teratorn is a unique herpesvirus that retains its viral genes intact and has acquired the endogenized lifestyle by hijacking the transposon system. However, it is unclear how this novel element evolved in the teleost lineage and whether fusion of two mobile elements is a general phenomenon in vertebrates. RESULTS: Here we performed a comprehensive genomic survey searching for Teratorn-like viruses in publicly available genome data and found that they are widely distributed in teleosts, forming a clade within Alloherpesviridae. Importantly, at least half of the identified Teratorn-like viruses contain piggyBac-like transposase genes, suggesting the generality of the transposon-herpesvirus fusion in teleosts. Phylogenetic tree topologies between the piggyBac-like transposase gene and herpesvirus-like genes are nearly identical, supporting the idea of a long-term evolutionary relationship between them. CONCLUSION: We propose that piggyBac-like elements and Teratorn-like viruses have co-existed for a long time, and that fusion of the two mobile genetic elements occurred frequently in teleosts.

Unlinking the methylome pattern from nucleotide sequence, revealed by large-scale in vivo genome engineering and methylome editing in medaka fish.

The heavily methylated vertebrate genomes are punctuated by stretches of poorly methylated DNA sequences that usually mark gene regulatory regions. It is known that the methylation state of these regions confers transcriptional control over their associated genes. Given its governance on the transcriptome, cellular functions and identity, genome-wide DNA methylation pattern is tightly regulated and evidently predefined. However, how is the methylation pattern determined in vivo remains enigmatic. Based on in silico and in vitro evidence, recent studies proposed that the regional hypomethylated state is primarily determined by local DNA sequence, e.g., high CpG density and presence of specific transcription factor binding sites. Nonetheless, the dependency of DNA methylation on nucleotide sequence has not been carefully validated in vertebrates in vivo. Herein, with the use of medaka (Oryzias latipes) as a model, the sequence dependency of DNA methylation was intensively tested in vivo. Our statistical modeling confirmed the strong statistical association between nucleotide sequence pattern and methylation state in the medaka genome. However, by manipulating the methylation state of a number of genomic sequences and reintegrating them into medaka embryos, we demonstrated that artificially conferred DNA methylation states were predominantly and robustly maintained in vivo, regardless of their sequences and endogenous states. This feature was also observed in the medaka transgene that had passed across generations. Thus, despite the observed statistical association, nucleotide sequence was unable to autonomously determine its own methylation state in medaka in vivo. Our results apparently argue against the notion of the governance on the DNA methylation by nucleotide sequence, but instead suggest the involvement of other epigenetic factors in defining and maintaining the DNA methylation landscape. Further investigation in other vertebrate models in vivo will be needed for the generalization of our observations made in medaka.

Centromere evolution and CpG methylation during vertebrate speciation.

Centromeres and large-scale structural variants evolve and contribute to genome diversity during vertebrate speciation. Here, we perform de novo long-read genome assembly of three inbred medaka strains that are derived from geographically isolated subpopulations and undergo speciation. Using single-molecule real-time (SMRT) sequencing, we obtain three chromosome-mapped genomes of length ~734, ~678, and ~744Mbp with a resource of twenty-two centromeric regions of length 20-345kbp. Centromeres are positionally conserved among the three strains and even between four pairs of chromosomes that were duplicated by the teleost-specific whole-genome duplication 320-350 million years ago. The centromeres do not all evolve at a similar pace; rather, centromeric monomers in non-acrocentric chromosomes evolve significantly faster than those in acrocentric chromosomes. Using methylation sensitive SMRT reads, we uncover centromeres are mostly hypermethylated but have hypomethylated sub-regions that acquire unique sequence compositions independently. These findings reveal the potential of non-acrocentric centromere evolution to contribute to speciation.

Hypomethylated domain-enriched DNA motifs prepattern the accessible nucleosome organization in teleosts.

BACKGROUND: Gene promoters in vertebrate genomes show distinct chromatin features such as stably positioned nucleosome array and DNA hypomethylation. The nucleosomes are known to have certain sequence preferences, and the prediction of nucleosome positioning from DNA sequence has been successful in some organisms such as yeast. However, at gene promoters where nucleosomes are much more stably positioned than in other regions, the sequence-based model has failed to work well, and sequence-independent mechanisms have been proposed. RESULTS: Using DNase I-seq in medaka embryos, we demonstrated that hypomethylated domains (HMDs) specifically possess accessible nucleosome organization with longer linkers, and we reassessed the DNA sequence preference for nucleosome positioning in these specific regions. Remarkably, we found with a supervised machine learning algorithm, k-mer SVM, that nucleosome positioning in HMDs is accurately predictable from DNA sequence alone. Specific short sequences (6-mers) that contribute to the prediction are specifically enriched in HMDs and distribute periodically with approximately 200-bp intervals which prepattern the position of accessible linkers. Surprisingly, the sequence preference of the nucleosome and linker in HMDs is opposite from that reported previously. Furthermore, the periodicity of specific motifs at hypomethylated promoters was conserved in zebrafish. CONCLUSION: This study reveals strong link between nucleosome positioning and DNA sequence at vertebrate promoters, and we propose hypomethylated DNA-specific regulation of nucleosome positioning.

Complete fusion of a transposon and herpesvirus created the Teratorn mobile element in medaka fish.

Mobile genetic elements (e.g., transposable elements and viruses) display significant diversity with various life cycles, but how novel elements emerge remains obscure. Here, we report a giant (180-kb long) transposon, Teratorn, originally identified in the genome of medaka, Oryzias latipes. Teratorn belongs to the piggyBac superfamily and retains the transposition activity. Remarkably, Teratorn is largely derived from a herpesvirus of the Alloherpesviridae family that could infect fish and amphibians. Genomic survey of Teratorn-like elements reveals that some of them exist as a fused form between piggyBac transposon and herpesvirus genome in teleosts, implying the generality of transposon-herpesvirus fusion. We propose that Teratorn was created by a unique fusion of DNA transposon and herpesvirus, leading to life cycle shift. Our study supports the idea that recombination is the key event in generation of novel mobile genetic elements. Teratorn is a large mobile genetic element originally identified in the small teleost fish medaka. Here, the authors show that Teratorn is derived from the fusion of a piggyBac superfamily DNA transposon and an alloherpesvirus and that it is widely found across teleost fish.